Impact ionization in InSb probed by THz-pump THz-probe spectroscopy 
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Picosecond carrier dynamics in indium antimonide (InSb) following excitation by below-bandgap 
broadband far infrared radiation were investigated at 200 K and 80 K. Using a novel THz- 
pump/THz-probe scheme with pump THz fields of 100 kV/cm and an intensity of 100 MW/cm^, we 
observed carrier heating and impact ionization. The number of carriers produced exceeds lO^^cm""^, 
corresponding to a change in carrier density AN/N of 700% at 80K. The onset of a well defined 
absorption peak at 1.2 THz is an indication of changes in LO and LA phonon populations due to 
cooling of the hot electrons. 
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Hot electron dynamics and carrier multiplication in 
semiconductors play important roles in the design of 
fast electronic devices as well as in the development of 
highly efficient solar cells [H. The elucidation of carrier 
dynamics on the ultrashort timescale is hence of great 
technological as well as fundamental interest. Indium 
antimonide (InSb) is a model system for the study of 
hot electron dynamics due to its low bandgap of 170 
mcV at room temperature and the fact that it has 
the highest electron mobility and saturation velocity of 
all known semiconductors. Its high mobility allows the 
fabrication of transistors with extremely high switching 
speed, but further technological applications are com- 
plicated by the low impact ionization threshold. Impact 
ionization by high electric fields is a well known phe- 
nomenon in InSb [J and has been studied mostly in 
the quasi DC limit of ujt <^ 1 where r is the charac- 
teristic electron momentum relaxation time, typically of 
the order of a few picoseconds. In this limit, the ex- 
tent of impact ionization is determined by the probabil- 
ity of an electron gaining enough energy from the driv- 
ing field to cross the ionization threshold. An equilib- 
rium between energy gain of the electrons due to the 
accelerating field and energy loss due to phonon scatter- 
ing is reached, leading to saturation in the drift velocity 
frf. Carriers with energies greater than the bandgap Sg 
can produce new electron-hole pairs through the impact 
ionization process, which also can be viewed as an in- 
verse Auger effect The newly generated carriers can 
be accelerated by the field and can induce impact ion- 
ization themselves. Carrier relaxation dynamics in the 
high-frequency limit (lut ^ 1) have been studied in InSb 
by sub picosecond time-resolved mid-infrared nonlinear 
spectroscopy using above-bandgap excitation Q . Below- 
bandgap excitation by CO2 lasers also has been used 
for time resolved measurements with nanosecond reso- 
lution Q and for intensity-dependent transmission mea- 
surements [1, Q . Below-bandgap excitation allows study 
of bulk hot-carrier dynamics without optical generation 
of holes or strong absorption and high carrier densities 



very close to the surface. 

In this letter, we report experimental observations of 
carrier generation in InSb at 80K and 200K due to im- 
pact ionization induced by below-bandgap IR radiation 
on the picosecond time scale, where ojt « 1. Near single- 
cycle pulses with fields strengths up to 100 kV/cm and 
a duration of 1 ps were used. The rise time of the THz 
pulses was less than the electron momentum relaxation 
time of T = 2.5 ps in InSb at 77 K. This leads to highly 
accelerated carriers that can cause carrier multiplication 
through impact ionization [loL 11|. THz radiation also 



can be used as a very sensitive probe to directly monitor 
free carrier behavior in semiconductors [l^ . The combi- 
nation of sub-bandgap direct excitation of doped semi- 
conductors and time resolved spectroscopy pro vides an 
excellent tool for observing carrier dynamics 13, l3, l3l ■ 
The experimental setup shown in Figure[l]was used for 
collinear THz-pump/THz-probe measurements. We gen- 
erated single-cycle THz pulses by optical rectification of 
a 800 nm, 5.5 mJ pulse from a TiiSapphirc laser with 100 
fs pulse duration at a repetition rate of 1 kHz. We tilted 
the pulse intensity front with a grating-lens combination 
to achieve noncoUinear velocity matching in lithium nio- 
bate 3.3, yielding THz pulses with energies up to 3 
The optical beam was split using a 10:90 
beamsplitter into two parts that were recombined un- 
der a small angle at the same spot on the grating. The 
10% part was modulated using an optical chopper and 
was used to generate the THz probe. The 90 % part 
was variably delayed and was used to generate the THz 
pump pulse. The single-cycle THz pulses were focused 
onto the sample using a 90-degree off-axis parabolic mir- 
ror pair with 190 and 75 mm focal lengths. The focused 
size was determined to be 1 mm using a razor blade 
scan. Another off-axis parabolic mirror pair with focal 
lengths of 100 and 190 mm was used to image the sam- 
ple plane onto the ZnTe detector crystal for electro-optic 
sampling of the THz field using balanced detection and 
a lock-in amplifier [l9| . In order to ensure the linearity 
of the detected signal and to eliminate THz pulse reflec- 
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FIG. 1: The experimental setup uses pump and probe THz 
pulses in a coUinear geometry, generated by tilted pulse front 
excitation in LiNbOa (LN) 



tions [201 , the ZnTc sampling crystal had an active layer 
of 0.1 mm and a total thickness of 1.1 mm. Selective 
chopping of the probe beam provided excellent suppres- 
sion of the pump pulse. Spectral analysis of our THz 
pump-probe results was conducted in the 0.2 to 1.6 THz 
range where the spectral amplitude was sufficiently high. 
A pair of wiregrid polarizers was used to attenuate the 
THz pulses for intensity-dependent studies. The samples 
were a n-type Te-doped and a nominally undoped InSb 
wafer, each 450 /xm thick, with carrier concentrations at 
77 K of 2.0 X lO^^cm-^ and 2.0 - 4.9 x lO^cm'^ respec- 
tively. The mobility as specified by the manufacturer was 
2.5 X lO^cm^/Vs. The THz fields were polarized paral- 
lel to the (100) axes of the crystals. We measured the 
THz fields E{t) that reached the ZnTe crystal with and 
without the sample in the beam path from which we cal- 
culated the effective absorption coefficient 



aeff = -3 In 
a 



(1) 



where d is the sample thickness, <max is the time win- 
dow of the measurement and T is a factor accounting 
for reflection losses at the sample surfaces. The quan- 
tity aeff is equivalent to the energy absorption coefficient 
averaged over our bandwidth. Figure [5] shows the time- 
resolved absorption traces for the doped sample at 80 
K and 200 K covering a probe delay range of up to 30 
ps. At both temperatures, the absorption increases after 
THz excitation and reaches a plateau after 30 ps, with 
a total increase in absorption of 80-90 cm~^. This rise 
is caused by the generation of new carriers through im- 
pact ionization. Unique to the measurement at 200 K 
is the initial dip in the absorption immediately after the 
THz excitation. The drop in absorption in this case is 
caused by a decrease in the mobility of the highly en- 
ergetic hot electrons as a result of both the strong non- 
parabolicity of the T valley in the conduction band of 
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FIG. 2: THz-pump/THz- probe time-resolved absorption 
data, spectrally averaged over the 0.2-1.6 THz range, from 
doped InSb [N^ = 2.0 x lO^^cm'^) at 80 K and 200 K. The 
deviation starting at 10 ps is due to re-heating caused by the 
THz pulse reflection within the sample. 



InSb |2l| and scattering of these hot electrons into side 
valleys [l; 



I22I . l23j . The carrier mobility /j, directly 
influences the absorption coefficient a which is typically 
conceptualized by the Drude model, parametrized by the 
plasma frequency Wp = (A'^e^/eoCoo'Ti*)^/^ and the mo- 
mentum scattering rate 7 in the form 



a 



nc{Lj'^ + 7^) ' 



(2) 



where is the high-frequency limit of the dielectric 
function and n is the refractive index. In the low- 
frequency limit {lo <C 7) the absorption a is directly 
proportional to the carrier mobility ^ — e/jm*. Due 
to the small bandgap of InSb, the intrinsic carrier con- 
centration at 200 K is much higher than that at 80 K. 
The effect of absorption saturation is thus much stronger 
at 200 K. After 4 ps, the additional absorption caused 
by the newly generated carriers exceeds the saturation 
effect, the magnitude of which is also diminishing as a 
result of the cooling of the hot carriers, thereby leading 
to the delayed rise observed in the overall absorption. 

Comparison between the equilibrium absorption of 
InSb at 80 K and the value measured 30 ps after intense 
THz excitation shows a eight-fold increase, indicating 
a similar increase in carrier concentration from 2x10^^ 
cm~'^ to 1.5 X lO^^cm"^. The same analysis cannot be 
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applied reliably at 200 K due to the very large intrinsic 
carrier absorption, which overwhelms the dynamic range 
of our spectrometer system. This effect is especially acute 
at low frequencies where the absorption is strongest, lead- 
ing to an apparent saturation of the frequency- averaged 
absorption shown in Figure [^b- Some signal at t < ap- 
pears in Fig. [2^ because of the nonlinear interaction in 
the LN crystal between the THz pump and probe fields 
and the optical pulses that generate them as well as non- 
linear interactions in the sample. 

We employ a simple system of rate equations, displayed 
below, to model the dynamics of impact ionization in the 
first 30 ps, assuming quadratic scaling Q of the impact 
ionization probability above the threshold energy eth ~ 

Eg. 

^ ^ C(e(i) - stk)^ ■ N{t) ■ eie{t) - (3) 

5 = - sth)^ ■ eth ■ e(e(t) - eth) ~ - (4) 

at 

In this model N{t) is the electron concentration, e(t) 
is the average carrier energy and 0(t) is the Heavi- 
side step function. We used a numerical value of C = 



7 X 10'^" J- 



obtained from Ref. [2J]. The energy 



relaxation time was assumed to be time and energy 
independent. 

A numerical solution, taking into account the reflec- 
tion at the sample interface, is shown in Fig. [31 The 
effect of absorption saturation due to carrier heating 
(isl . [l^ was accounted for approximately by assuming 
a^jj oc iV(i)[eo — e(i)]/eo where eg is the average car- 
rier energy immediately after excitation. From this fit 
we obtain a phenomenological relaxation time time 
of 7 ps, much longer than the 1.3-2 ps calculated for 
the energy relaxation time in the DC limit by Kobayashi 
[2^. The faster value cannot describe our experimental 
conditions since it would preclude impact ionization and 
carrier cooling dynamics that we observe at slower time 
scales. In order to elucidate the lattice dynamics further, 
we performed a series of intensity-dependent pump-probe 
measurements with a fixed probe delay. Figure |4] shows 
the absorption spectra obtained with different pump in- 
tensities at a probe delay of 35 ps for the doped and un- 
doped samples at 80K. At frequencies below 0.6 THz, we 
observe the expected Drude-like contribution from free 
carrier absorption which is more pronounced at higher 
pump fluence. In addition, we observe a distinct ab- 
sorption peak at 1.2 THz in the undoped sample and a 
weak feature that indicates a similar peak in the doped 
sample. The amplitude of the peak is highly intensity de- 
pendent, and appears to approach its asymptotic value 
just above 50% of the maximum intensity. The behavior 
of this peak suggests that its origin is lattice vibrational 
rather than electronic. Polar optical phonon scattering 
is well known as the dominant energy loss mechanism for 
hot electrons in InSb pjj . The main channel of energy 
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FIG. 3; Experimental data for the undoped sample (squares) 
and simulation results for the carrier concentration (dashed) 
and absorption (solid line) based on quadratic scaling of im- 
pact ionization rate with carrier energy. Parameters used 
were: N(t = 0) = 5 x 10^^ cm~^, energy relaxation time 
Te = 7 ps and eo = e(< = 0) = 1.3 eV. 



loss is through the emission of LO phonons with a fre- 
quency of 5.94 THz. These phonons decay into acoustic 
modes through anharmonic terms of the crystal potential 
and through the second-order electric moment of the lat- 
A series of sum- and difference phonon peaks 
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tice 

between 1 to 10 THz has been observed and assigned 
in the far- infrared spectrum of InSb [2^. At very low 
THz fields, produced by a photoconductive antenna, we 
also were able to observe some of these weak absorption 
peaks in the undoped sample. The assignments reported 
in [2^ indicate a 1.2 THz difference frequency between 
the LO and LA mode at the zone boundary. The dras- 
tic change in the absorption coefficient of the difference 
phonon peak is the result of large changes in phonon pop- 
ulations due to energy transfer from the hot electrons 
generated by the THz pump pulse. Monte-Carlo simula- 
tions [29! have shown that substantial phonon population 
changes can occur even at comparatively low DC fields 
on picosecond timescales. The temporal evolution of the 
electron-lattice interaction can be studied by separating 
the spectrally resolved pump probe data into two differ- 
ent frequency bands 0.2-0.6 THz and 1.0-1.2 THz. This 
is illustrated in Figure [5] where data from time-resolved 
measurements at full THz intensity and at half the full 
intensity are shown for the doped sample. For the low- 
frequency band, the rise in absorption is almost identical 
at the two intensities, while the contribution from the ab- 
sorption between 1.0 to 1.2 THz is reduced considerably 
when the THz intensity is halved. The higher frequency 
portion of the absorption spectrum is clearly more sensi- 
tive to the intensity of the THz pump pulse, as illustrated 
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FIG. 4: THz absorption spectra at various pump intensities 
measured at a probe delay time of 35 ps for the doped sample 
(a) and undoped sample (b) at 80K 
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FIG. 5: Time-resolved spectrally averaged absorption for the 
doped sample at 80 K with 100% pump intensity (black 
squares) and 50% pump intensity (red triangles), (a) the av- 
erage absorption between 0.2 and 0.6 THz, (b) the average 
absorption between 1.0-1.2 THz 



in Figure[4jD. In addition, at the lower intensity there is a 
delay of approximately 10 ps before the rise of the absorp- 
tion signal. This delay was reproduced qualitatively for 
all measurements at intermediate THz pump levels. We 
do not have a complete understanding of the intensity- 
dependence of the delay, but we note that the delay is 
on the order of the LO phonon decay time in InSb 27| , 
suggesting that it arises from electronic and LO phonon 
relaxation processes that populate the LA phonons. The 
large fluctuations observed at 12 and 24 ps are due to 
THz pump pulse reflections in the sample that are over- 
lapped in the EO crystal with the THz probe pulse. 

The newly developed THz pump/THz probe technique 
permits sensitive monitoring of carrier dynamics in semi- 
conductors on the picosecond timescale. We observed 
the dynamics of impact ionization and carrier generation, 
up to a seven-fold increase over the equilibrium carrier- 
generation, following intense THz excitation of InSb. Our 
ability to spectrally analyze the time-resolved signal al- 
lowed us to detect distinct features that we attribute de- 
cay of electronic energy into LO and LA modes. Monte- 
Carlo simulations are needed for a more in-depth un- 
derstanding of the interplay between hot electrons and 
lattice, taking into account effects like impact ioniza- 
tion, intervalley- and polar optical phonon scattering and 
changes in phonon population. Additional effects of the 
THz fields themselves, including the possibility of THz- 
induced band-to-band tunneling [sot that could produce 
new carriers directly, also warrants further analysis. 
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